Helium line formation and abundance during a C-class flare 
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1 — '; ABSTRACT 

7-H ■ 

qq ■ During a coordinated campaign which took place in May 2001, a C-class flare 

was observed both with SOHO instruments and with the Dunn Solar Telescope 
of the National Solar Observatory at Sacramento Peak. In two previous papers 
we have described the observations and discussed some dynamical aspects of the 
■ earlier phases of the flare, as well as the helium line formation in the active region 

prior to the event. Here we extend the analysis of the helium line formation to 
the later phases of the flare in two different locations of the flaring area. We 

^ ■ have devised a new technique, exploiting all available information from various 

o5 ■ SOHO instruments, to determine the spectral distribution of the photoionizing 

EUV radiation produced by the corona overlying the two target regions. In order 
to find semiempirical models matching all of our observables, we analyzed the 
effect on the calculated helium spectrum both of A He (the He abundance) and of 
the uncertainties in the incident EUV radiation (level and spectral distribution). 
We found that the abundance has in most cases (but not in all) a larger effect 
than the coronal back-radiation. The result of our analysis is that, considering 
the error of the measured lines, and adopting our best estimate for the coronal 
EUV illumination, the value An e = 0.075 ± 0.010 in the chromosphere (for T > 
6300 K) and transition region yields reasonably good matches for all the observed 
lines. This value is marginally consistent with the most commonly accepted 
photospheric value: An e = 0.085. 
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Introduction 



A key ingredient to understand the physics of solar and stellar plasmas is their chem- 
ical composition. An astrophysical plasma is subject to c hemical frac tionation processes 
of diverse nature (gravity, thermal diffusion, etc., see e. g. iDrakel 120031 ) that can produce 
abundance "anomalies" between regions of different temperature and density. For the solar 
corona these anomalies are generally described in terms of the so-called FIP (First Ionization 
Potential) effect. Not only the details, but the very general physical framework behind this 
effect is still debated. However, one of the few aspects over whi ch there is appar ently a 
consensus, is that the FIP effect should arise in the chromosphere (lGeisslll982l . Il998l ). 



Helium is one of the few elements that exhibits strong lines forming in the chromosphere, 
and thus, in principle at least, its abundance ([He/H] = logA He (where A He = N He /N H , the 
ratio of number densities of He and H) could be estimated in that region. It is also the 
element with the highest FIP (24 eV). Furthermore, it is the second most abundant element 
in the Sun (and in the universe, for that matter), and therefore must be included in any 
theoretical model of fractionation processes in the solar atmosphere. Moreover, the coronal 
[He/Hi strongly depends on the availability of helium in the upper chromosphere (e. g.: 
Killie et al.l[2005h . 



Unfortunately, analyzing helium lines in the solar atmosphere is a particularly challeng- 
ing task, more so than other "typical" chromospheric lines such as Ha o r the Ca II H fc K 
doub l et. The main reasons ar e described in several papers (for example: lAndretta fc Jones 
19971 ; iPietarila &: Judge! 120041 ) . Among the problems that can be mentioned: all of its lines 
need detailed radiative transfer calculations; strong interlocking between n = 2 and n = 3 
singlet and triplet levels makes the simultaneous treatment of EUV and optical lines almost 
unavoidable; the high excitation energies of all helium lines (a consequence of its high FIP 
value) makes those lines very sensitive to non-equilibrium effects; the ionization rates of neu- 
tral and ionized helium can be significantly affected by EUV radiation via the so-called "Pho- 
toionization - Recombination" (P-R) mechanism (e. g.: IZirinlll975l ; lAndretta et al.ll2003l ). and 
therefore not only coronal plasmas, but even the strong He II Ly-a at 30.4 nm can be signif- 
icant photoionization sources for neutral helium (another example of the strong interlocking 
effects in the helium atomic system). 

All these problems imply, among other things, that a realistic study of helium lines in 
the solar chromosphere requires quite sophisticated modeling, as well as simultaneous mea- 
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surements in both the EUV and the optical (visible and IR) range. To this aim we planned an 
observing campaign (SOHO Joint Observing Programme 139) coordinated between ground 
based and SOHO instruments to obtain simultaneous spectroheliograms of the same area in 
several spectral lines, including four He lines (He I 58.76, He I 1083.0, He I 58.4 and He II 
30.4 nm), that sample the solar atmosphere from the chromosphere to the transition region. 
The EUV radiation in the range A < 50 nm and in the range 26 < A < 34 nm has also been 
measured at the same time. 



In a previous paper (jMauas. Andretta. Falchi. et al.ll2005l . hereafter Paper II), we carried 
out an analysis of the helium line formation in a specific location of the observed active region 
(NOAA AR 9468, cos6>=0.99), prior to a small two-ribbon flare (GOES class CI) developed in 
this region on 2001, May 26, around 16:00 UT. We found that the incident coronal radiation 
ha s a limited effect on the UV He lines (in partial agreement with the proxy-based analysis 
of iFredvik fc Maltbyill999l . of the He I 58.4 line), while being fundamental for the optical 
lines. For these lines, we also found that photons of shorter wavelengths are more effective 
at increasing their line depth. Incidentally, this latter finding implies that for the calculation 
of the He line profiles it is important, at least in principle, to correctly estimate not only the 
total number of ionizing photons, but also their spectral distribution. But the main result 
of Paper II was that the photospheric value of [He/H] is compatible with our observations 
within an uncertainty of about a factor two and that such uncertainty could be reduced in 
a higher density regime. 

On the basis of these results we decided to extend our analysis to the data taken during 
the flare. In this paper we present results of such analysis and study the formation of He 
lines after the main evaporation phase, when the He lines presumably originate in a high 
density plasma and the incident coronal radiation, larger than in the pre-flare phase, can 
be estimated with a lower uncertainty. We construct semiempirical models of the flaring 
atmosphere to match the observed line profiles from the chromosphere to the transition 
region, taking into account the total EUV irradiance in the range A < 50 nm. We also 
analyze the effect of [He/H] changes, with the purpose of attaining an empirical estimate of 
that parameter. 

To this regard, we mention that techniques based on 7-ray line measurements produced 
in flares offer an alternat i ve ayenue for estimating the chromospheric He abundance (e. g. 
Mandzhavidze et al.lll997l . ll999l ; IShare fc Murphylll998l ). Results from these techniques seem 
generally to indicate an enhancement either of the accelerated a/proton ratio or of the 
"ambient" He abundance. Attempts at resolving this a mbiguity seem however to hint that 
^He is consistent with its accepted photospheric values (jMurphyl 120071 ). 



For completeness, we should also mention recent EUV spectroscopic measurements of 



-4- 



A\\ e , during a flare, m ade by analyzing the Balmer-7 of He II at 108.5 nm in off-limb spectra 



( IFeldman et al.ll2005l ). which yielded Au e ~ 0.122 ± 0.024, marginally higher than the pho- 
tospheric value. This result, however, refers to off-limb post-flare loops and thus cannot be 
directly compared with our on-disk, chromospheric measurements. 

In Sec. [2] we briefly describe the set of observations, focussing our attention on two 
specific locations of the flaring region. In Sec. [3] we describe the procedures employed to 
estimate the incident coronal EUV radiation illuminating the regions under examination and 
in Sec. H] the chromospheric modeling for those two locations. Finally, in Sec. [5] we explore 
how the helium spectrum is affected by the value of the chromospheric He abundance and 
by the uncertainties in the estimates of coronal radiation. 



The Observations 



A detailed description of the observing program is given in lTeriaca. Falchi. et al.l (120031 ) 
(hereafter Paper I). We recall here that spectroheliograms were acquired with the Horizontal 
Spectrograph at the Dunn Solar Telescope (DST) of the National Solar Observatory at Sa- 
cramento Peak in the chromospheric lines Ca II K, Ha, and Na I D as well as in the He I 
lines at 587.6 (D 3 ) and 1083 nm. The full field of view (FOV) of the DST, 170" x 170", was 
covered in about 5 minutes, with a sampling step of 2". Correcting for offsets among the 
different detectors resulted in a final useful FOV of 160" x 140" with an effective resolution 
of 2". 

During the same period, spectroheliograms of the active region were obtained in selected 
spectral lines wit h the Normal Inciden ce Spectrometer (NIS) of the Coronal Diagnostic Spec- 
trometer (CDS) ( Harrison et al.lfl.995l ) aboard SOHO. The selected lines are: Fe XVI 36.1 nm 
in the NIS band #1; He I 58.4 nnu He II 30.4 nm (2 nd order), O V 62.9 nm, and the blend 
Fe XII 59.26 + Fe XIX 59.22 nnfl in the NIS band #2. The 4"-wide slit was stepped 6" 
covering a 148" wide area in ~ 5.5 minutes. The final useful FOV was 148" x 138" with an 
effective pixel size of 6" x 3". 4. Ground-based and CDS data were aligned using SOHO/MDI 
images as a reference. We estimate the error around few arcseconds. CDS and ground based 
spectra are simultaneous within 2 minutes. 



Another SOHO instrument, the Extreme ultraviolet Imaging Telescope (EIT; iDelaboudiniere et al 



1 In the case of the Fe XII 59.26+Fe XIX 59.22 nm blend, an analysis of the mean line wavelengths shows 
(see also Paper I) that the former contribution is small or even negligible in the areas interested by the flare, 
but, conversely, is dominant in the surrounding, unperturbed regions. 
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19951 ) provided synoptic series of full disk images centered on 17.1, 19.5, 28.4 and 30.4 nm 
around 13:00 and 19:00 UT, plus a series of full disk images centered on 19.5 nm acquired 
with a 12 minutes cadence during the routine CME watch program. Fig. [T] shows a few 
images extracted from the EIT 19.5 sequence, showing the development of the flare. 



16:00:10 UT 16:12:10 UT 16:24:11 UT 16:36:11 UT 



16:48:11 UT 
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-100 -50 
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Fig. 1. — The evolution of the NOAA AR 9468 in the first hour of the event, as seen in the 
sequence of SOHO/EIT 19.5 nm images (see Fig. [2] for a comparison with XUV irradiance 
and Ha radiance light curves). The contours represent Ha isophotes relative to the DST 
raster scan nearest in time. The positions of the two regions studied in detail in this paper 
(regions A and D) are marked with crosses. 

Finally, the EUV flux during the flare (see Figure El) was ob served by the Solar EUV 



Monitor (SEM) instrument aboard SOHO flHovestadt et al.lll995f ). The CELIAS/SEM in- 



strument provides calibrated total photon counts in the range A < 50 nm (zero-th order, 
SEM ), and in the range 26 < A < 34 nm (first order, SEMi), at 1 AU. 

The flare dynamics has been discussed in Paper I, where it has been shown that a large 
area of the flaring region is affected by chromospheric evaporation in a short time interval 
(16:02:30 - 16:05:30 UT) during the first impulsive phase. In the present paper we analyze 
the He line profiles after 16:13 UT, when the chromospheric evaporation is almost finished, 
in two distinct locations, which we term regions 'A" and "D" (size of 4" x 2" each). These 
two regions are marked in Fig. [1] 

From what we obtained in Paper I, region A did not undergo any evaporation process, 
at least within our limit of detectability, while region D is in the area where we measured 
downflows at chromospheric levels and upflows at Transition Region (TR) and coronal levels. 
Moreover, while in region A the chromospheric radiance is almost constant after 16:13 UT, 
in region D it approximately follows the EUV light curve (see Figure [2]) with a secondary 
maximum at about 16:20 UT, showing that also this small flare is a multiple event. A similar 
behaviour is also detected for the TR and coronal lines with temperatures up to 10 6 K. 

We use region A as a control region where only small changes are detected in our 
observed features. We construct (Sec. H]) semiempirical models of the flaring atmosphere for 
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this region at two times when the SEM EUV irradiance is nearly the same: 16:14 and 16:33 
UT (dashed lines in the SOHO/SEM light curve in the upper panel of Fig. [2]). 

For region D the models are constructed at 16:15, 16:20, and 16:35 UT (dotted lines). 
These times will be also referred to as t ly t 2 , and t 3 , respectively. Time t 2 corresponds to 
a peak of the total SEM irradiance, while at t\ and t% the irradiances are almost equal. 
Therefore, the choice of these regions and times allows us to check whether and how the He I 
lines are influenced by the EUV coronal back-radiation, and/or by the atmosphere where 
they form. This point will be discussed more in detail in Sec. 13.11 

The two sets of optical line profiles are shown in Fig. [3], with the exception of the Na I 
doublet that did not change during the flare in the two considered regions. For region A, 
since the profiles remain practically unchanged, only the profiles observed at 16:14 are shown 
and compared with profiles taken in a nearby, relatively quiescent region, prior to the flare, 
which have been used to calibrate the spectra in absolute units. 

The noise of individual spectra is sufficiently small to be negligible compared to the pixel- 
to-pixel variability within the target regions which, therefore, determines the uncertainty on 
the mean profiles. The uncertainties as function of wavelength on the mean profiles in 
region A at 16:14 UT and in region D at 16:20 UT are indicated in Fig. [3] as a grey band. 
In practice, those bands are obtained, at each wavelength, from the r.m.s. of the iV profiles 
within the target regions, divided by y/~N. 

It is important to realize that for the two relatively weak He I lines, the uncertainty in 
the central line depth (the quantity which carries most information about the chromosphere) 
is also due to the variability of the underlying photosphere. This fact is most evident in the 
continuum bands of the D 3 line at 58.76 nm. By examining the variability of the profiles 
normalized to the continuum, we estimate that for the latter line the photospheric variability 
contributes by about 50% to the variability of the central line intensitjIE In any case, the 
uncertainties shown in Fig. [3] for the two optical He I lines (of the order of 3%), are to be 
taken as upper limits for the uncertainties due to the variability of the chromosphere only. 

In the case of the EUV He I lines, the width of the instrumental response of the CDS 
spectrograph prevents a detailed comparison of computed and observed line profiles. We thus 
considered only the radiances. For these lines, the main source of error is the uncertainty in 



2 We adopt here, as commonly done in Astrophysics, the term "specific intensity" or simply "intensity" in 
place of the more formally correct "spectral radiance" (SI units: W sr _1 m~ 2 nm -1 ) and indicate it with the 
symbol I\; with "radiance" it is then meant an intensity integrated over a wavelength band or a line profile: 
f I\ dA (SI units: W sr _1 m~ 2 ). Finally, "irradiance" is the spectral and angular integral of an intensity: 
fix dA dfl (SI units: W m" 2 ). 
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the absolute radiometric calibration, estimated to be of the order of 30%. 

3. Coronal back-radiation 

In Paper II the average EUV illumination over the target area in active region AR 9468, 
integrated in the A < 50 nm range, was computed using the spatial information provided 
by the EIT synoptic images to estimate the contribution to the full-disk SEM 0.1-50 nm 
irradiance due to the target region. This value, J re f= 1.2 ±0.2 x 10 15 photons s _1 cm -2 sr _1 , 
is considered here as the reference value. We recall that the emission from the He II lines 
has been subtracted from 7 re f because these lines are taken into account in a self-consistent 
way in the model computations. 

Following the same approach, and using the EIT 19.5 images as proxies, we estimated 
that, in the hour before the flare, the EUV radiation in regions A and D was, in average, 
about 2.2 and 1.6 times higher than 7 re f. Since there are variations of the order of 30% in 
a 3 x 3 pixel box (« 16" x 16") around the center of each region, the standard deviation of 
the mean value in the same box would be 3 times smaller. We however prefer to use the 
3-er value on the mean value quoted above, as it is large enough to encompass other possible 
sources of error (for instance those introduced by the use of EIT 19.5 images as proxy for 
the full 1-50 nm range). Thus, we estimated the EUV illumination over regions A and D to 
be, respectively, about 2.6 x 10 15 and 1.9 x 10 15 photons s~ x cm -2 sr _1 (Tabled]). 

The implicit assumption in such an approach, as explained in Paper II, is that the 
spectral distribution of the EUV radiation over the active region does not change significantly. 
For the observations taken during the flare, however, this assumption may not be valid. Thus, 
we adopted an altogether different approach for estimating the EUV illumination at the times 
of our observations. 

From the EIT sequence of images we can estimate that the only event that can cause 
the increase in EUV irradiance shown in the topmost panel of Fig. [2J is the C-class flare we 
observed. From the increase in SEM irradiance, Af, and from the knowledge of the angular 
area affected by the flare, AQ, we can then estimate the increase in radiance in the flaring 
region, AI = Af/AQ. The total radiance in the region is then obtained by adding to AI 
the estimates of the pre-flare average radiances given above. 

For a first, rough estimate, we can start from a visual determination, from Fig. (2J of the 
excess irradiance in the SEMo measurements due to that event: Af « 2 to 4 x 10 8 photons 
s _1 cm -2 in the first 60 minutes of the event. On the other hand, from both EIT and CDS 
observations, we estimate flare areas to be of the order of AQ « 1500 to 2000 (") 2 in the 
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same time interval . Consequently, we expect radiances in the flaring regions to increase by 
something of the order of AI « 6 to 8 x 10 15 photons s -1 cm" 2 sr" 1 . 

The inaccuracy of this estimate is mainly due to the approximate definition of the flaring 
area: it can vary by as much as a factor 2 or 3, depending on which spectral signature is 
considered. For a more precise determination of the excess irradiances in regions A and D, 
we refined the method by using the information given by the variation in time (light curves) 
above the pre-flare mean value of both the SEM irradiance and of the observed spectral 
signatures in our data set. More in detail, we proceeded as follows: 

Step #1: We created a set of reference images, both from CDS and from EIT, by 
averaging the (co- aligned) images in the 40 minutes preceding the onset of the event. The 
variability of these images around the average is of the order of 0.1-0.2% (standard deviation), 
or even less. 

Step #2: We obtained maps of "residual" or "excess" radiance in each spectral feature 
during the flare, AI^(th), by subtracting the reference pre-flare images (taking into account 
solar rotation) for each time th, and for each pixel (i, j). Multiplying by the pixel angular size, 
dco k , we then obtain the contribution of each line and of each pixel to the EUV irradiance, 
A.fiJth) = AIh(th) du k . The result is shown in Fig. 0] for CDS images only. A similar 
procedure has been applied to the EIT sequence as well. More details on the analysis of the 
individual CDS rasters can be found in Paper I. Here we only mention that the rasters were 
taken by moving the slit from W to E — right to left, in the images of Fig. HJ where the 
beginning of each raster is indicated at the top of the corresponding column. 

Step #3: We computed the light curves of the excess irradiance for each spectral fea- 
ture by summing the contributions Affj(th). Each curve represents the contribution of the 
considered feature to the EUV excess irradiance during the flare, under the assumption that 
no other significant event is occurring on the solar surface outside the CDS FOV. We verified 
the correcteness of that assumption by inspecting the full-disk EIT 19.5 images during the 
event. The resulting light curves are shown in the top panel of Fig. [51 The set of CDS and 
EIT curves cover a broad range of temperatures (from logT < 5.0 for He I and He II, to 
logT ~ 6.9 for Fe XIX), and can thus be regarded as a good set of proxies for the temporal 
and spatial evolution of the various temperature components of the flaring plasma. 

Step #4: We found a linear combination of those light curves matching reasonably well 
the evolution of SEM measurements during the flare. More in detail, we modeled the total 
excess irradiance at each time th, Af(th), with a linear combination of the excess irradiance 
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of each line k, in the form: 

Af(t h ) ^a + ai t h + X) 6fc £ A 4(^) • W 

k ij 

The term a^th, that accounts for possible overall changes in the active region structure 
(hence: of the overall EUV emission) in the course of the event, is actually rather small: To 
reduce the number of free parameters, we set this term to zero. 

The two resonance helium lines, He I 58.4 and He II 30.4, exhibit an almost identical 
behaviour: only the latter has been taken into account in the linear regression. 

The comparison of the thus modelled EUV SEM irradiance (grey line in Fig. [5]) with the 
actual measured signal (smoothed to account for the 5-minutes duration of the CDS rasters) 
is shown in the bottom panel of Fig. [5J for the first ^120 minutes of the event. We stress 
that bk Ylij Af]* is not the contribution of the individual line k to the EUV irradiance, but 
is instead the contribution of the group of lines of similar temperatures, for which line k is 
a proxy. 

Beyond the mere application of the mathematical linear regression procedure, we think 
useful to remark here how the different parts of the irradiance SEM light curve are due to 
different temperature contributions of the flare plasma. More specifically, the first peak, at 
the time resolution of the CDS rasters, is modelled essentially by the cooler temperature lines 
(proxies: He II and O V), while the second peak (at ~16:20 UT) is due mostly to the hottest 
flaring plasma (proxy: Fe XIX), with a tail modelled by progressively cooler components 
(Fe XVI and then EIT 19.5). These finding, based on the modelling of the SEMo light curve, 
are qualitatively consistent with the characteristics of the other irradiance curves shown in 
Fig. [2J the second peak is indeed absent in SEM l5 while is more pronounced in the GOES 
signal. 

Step #5: The coefficients of the linear combination, bk, were then applied to the corre- 
sponding residual images on a pixel-by-pixel basis, producing a map of the contribution to 
the total excess EUV irradiance in the SEM spectral range (A < 50 nm) as function of time 
during the flare. The contribution due to transition region and coronal lines (henceforth, for 
brevity: "the coronal component"), needed in the model calculations, is obtained by sub- 
tracting the contribution of the proxy He II 30.4 nm, that will be computed self-consistently 
in the models. Thus, the contribution of each pixel to the EUV coronal irradiance, divided 
by the pixel's angular size, yields an estimate of the spatial variation of coronal radiance 
enhancement, in the band A < 50 nm, over the flare area. The result of such a procedure is 
shown in Fig. [6l 

If we make the reasonable assumption that the EUV emission during the flare is optically 
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thin, the maps of EUV excess radiance of Fig. [6] are also maps of EUV back-radiation 
illuminating the chromosphere. In Sec. 13.11 we will discuss more in detail this assumption. 

Step #6: Finally, by adding the pre-flare, average value of radiance estimated in the two 
target regions, A and D, we obtained an estimate of the value of EUV radiance of coronal 
origin in the band A < 50 nm, J cor , as function of time over the two target regions, and in 
particular at the times of the optical observations of Fig. [3] (Tabled]). 



Region 


Time 


Coronal radiance 




(UT) 


Excess 


Total 


Total// re f 


A 


pre-flare 




2.6 


2.2 




16:14 


0.38 


3.0 


2.5 




16:33 


1.73 


4.4 


3.6 


D 


pre-flare 




1.9 


1.6 




16:16 


5.48 


7.3 


6.2 




16:20 


7.38 


9.2 


7.7 




16:35 


5.25 


7.1 


6.0 



Table 1: Coronal radiances for regions A and D, in units of 10 15 photons s -1 cm -2 sr _1 . The 
value J re f = 1.2 x 10 15 is the reference (average pre-flare) value for the target area, estimated 
in Paper II. 

The (statistical) uncertainties in the excess number of photons are of the order of ~ 
0.4 x 10 15 (region A) and ^2x 10 15 (region D) photons s _1 cm -2 sr _1 , estimated as the 3-cr 
variation of the illumination around the target regions. 

We have also explored the systematic variations induced by changing the most relevant 
parameters that could affect the result of the procedure described above (for example: using 
a constant instead of a linear background in the regression calculations of Step #4); the 
resulting overall variations are smaller than the statistical uncertainty at the 3-cr level quoted 
above. 

Together with the uncertainty in the pre-flare radiance enhancements estimated earlier 
in this section (« 30%), and with the uncertainty of 7 re f estimated in Paper II (0.2 x 10 15 , or 
about 17%), we thus obtain an overall uncertainty in the total coronal radiance of the order 
of 25% to 35% in both regions A and D. 
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3.1. The spectral and angular distribution of the coronal back-radiation 



In Paper II we estimated the spectral distribution of the coronal photoionizing EUV 
radiation by taking a reference spectral distribution, and multiplying that distribution by 
a factor (close to unity) such that the radiance in the range A < 50 nm equals J re f = 
1.2 x 10 15 photons s" 1 cm" 2 sr _1 . 

In that paper we adopted as our reference spectral distribution the one for a "typical" 
active region comp uted with the spectral code and atomic database CHIANTI version 4 
( lYoung et al.l 120031 ) . The CHIANTI calculations require a description of the density and 
temperature distribution in the atmosphere via the Differential Emission Measure (DEM), 
defined as NuN e (dT / dh) _1 , where h is a coordinate along the line of sight; we adopted the 
table in file active regi on. dem provided with the package, which is based on data from 
Vernazza fe Reeved (119781 ). Furthermore, we adopted a "hyb rid" set of abundances ( inter- 
mediate between photospheric and coronal elemental mixtures: iFludra fe Schmelzlll999l ). and 
a constant pressure of P e /k = 3 x 10 15 cm -3 K. We verified that the use of the more re- 



cent version 5 of the same code (ILandi et al.l 120061 ) did not alter significantly that reference 
spectral distribution (Fig. 1 of Paper II). 

In the case of the flaring atmosphere we are considering, however, it is possible that 
significant changes in the spectral distribution may occurr. We tested this possibility by 
considering the case where departures from the "standard" spectral distribution (the pre- 
flare AR EUV spectrum, in fact) could be more pronounced: the emission above region D 
at 16:20 UT, where we observe a peak of total EUV emission (Tabled]), as well as a peak in 
the hot Fe XIX line (Fig. EJ. 

To do so, we found a DEM producing the measured total number of photons in the range 
A < 50 nm (9.2 x 10 15 photons s _1 cm~ 2 sr _1 ), and giving at the same time a reasonable match 
to the transition-region and coronal lines observed by CDS: O V 62.9 nm, Fe XVI 36.1 nm, 
and the blend Fe XII 59.26+Fe XIX 59.22 nm. The optically thick helium lines have not 
been considered in this analysis, but are nonetheless shown in the plots summarizing the 
results of this test (Fig. [7]). 

In the calculations, we assumed a higher pressure, P e /k = 3 x 10 16 cm -3 K, than in the 
pre-flare atmosphere. This value of the pressure corresponds, for instance, to N e « 10 11 and 
10 10 at logT = 5.5 and 6.5 respectivel y. These densities a re generally consistent with flare 



densities found in the lite rature (e. g. : Dere fe Cookl Il979l ) . although lower densities have 



sometimes been reported (ICook et al.l I1994J ) . Aside of the overall scale factorized by the 
DEM, the specific value for P e /k affects only the relative radiances of the density-dependent 
lines, and thus it is of little importance for the purpose of estimating the global flare spectral 
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distribution for A < 50 nm. 

Such a DEM is compared in the left-hand panel of Fig. [7] with the DEM used to compute 
the pre-flare, reference spectrum. The corresponding observed line radiances are shown i n 



that panel multiplied by the factor DEM(T e fj)/J ca i c (as in, e. g., iDel Zanna et al.ll2002l ); 
here T e g is the temperature of peak line emission. The error bars in the case of the data 
taken during the flare reflect the statistical variability in the 9 pixels around and above the 
target (D) region. Thus, we do not show the theoretical uncertainties or the calibration 
uncertainties both on the absolute intensities and on the CDS NIS1/NIS2 relative radiances 
(i. e.: the radiance of the Fe XVI line relative to all the others). 

The DEM-derived pre-flare and the flare spectra are compared in the right-hand panel 
of the same figure (dark and light grey histograms, respectively). We also show the pre-flare 
"standard" spectrum multiplied by the factor required to obtain the same total number of 
photons below 50 nm as in the flare spectrum (dashed histogram). The largest differences 
between these two spectra occur in the region below 15 nm, where the number of photons 
in the modelled flare spectrum is about a factor 4 higher than in the AR scaled spectrum. 
Conversely, the number of photons in the range 15-50.4 nm is about 20% smaller than in 
the AR scaled spectrum. 

On the other hand, for both distributions the bulk of the photoionizing photons is 
in the range 15 < A < 50 nm (about 3/4 of the total). In fact, this is the main reason 
why, as discussed in Paper II, the shorter wavelength part of the spectral distribution has 
relatively little effect on the photoionization of helium: We found that only the He I 1083 
line shows small changes due to the different penetration depths in the chromosphere of the 
EUV photons of different energies. We will further discuss this point in Sec. 15.21 

In conclusion, we expect, and indeed find (see Sec. I5.2p . that for our calculations it is not 
necessary to take into account the details of the distribution of photon energies: after all, we 
are studying a quite small flare. We thus simply scaled the "standard" spectral distribution 
by the appropriate factors listed in the last column of Table [H 

Concerning the problem of the angular distribution of the EUV photoionizing photons, 
we have already showed in Paper II that, despite their relatively featureless spatial distribu- 
tion across the target AR prior to the flare, at the chromospheric depths the details of the 
angular distribution are relatively unimportant compared to the overall spectral distribution 
and, most importantly, to the total number of impinging photons. 

In the case considered in this paper, on the other hand, during the flare the EUV 
emission is clearly concentrated in smaller areas (Figs. H]and[6]). It is therefore reasonable 
to assume that the main sources of ionizing radiation are concentrated directly above, or 
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very near the target regions. We verified this assumption using the control region (region 
A): during the time span we are considering, the helium lines in that region hardly change 
(Fig. [3]), despite the large variations of the EUV radiation around (above) region D (see 
Table [I]), which is roughly 40 Mm away. So, it is clear that the flaring chromosphere in 
region A, and a fortiori in region D, is sensitive only to the EUV radiation coming from 
regions significantly closer than 40 Mm pa 50", or, in other terms, is sensitive to the radiation 
coming at most from a few EIT or CDS pixels away. 



4. Chromospheric modeling 



Once the coronal back-radiation has been determined, semi-empirical models were con- 
structed to_jrj^tch_the_^served profiles. The modeling was done using the program PAN- 
DORA (lAvrett fe Loeserlll984l ). Given a T vs. h distribution, we solved the non-LTE radia- 
tive transfer and the statistical and hydrostatic equilibrium equations, and self-consistently 
computed non-LTE populations for 10 levels of H, 29 of He I, 15 of Fe I, 9 of C I, 8 of Si I, 
Ca I and Na I, 6 of Al I and He II, and 7 of Mg I. In addition, we computed 6 le yels of Mg II, 



and 5 of Ca II. The atomic models we used for H and Ca II are described in Mauas et al. 



(119971 ) and lFalchi fe Mauaa (11998I ). The models we used for He I and He II are described in 
Paper II. More details on the modeling , and on the different assumptions and their validity 
can be found in iFalchi fe Mauad (119981 ) and in Paper II. The adopted microturbulence, in 
particular, is the same as in the latter work. 

The modeling was done in a plane-parallel atmosphere. This is certainly an adequate 
approximation if we consider the much larger horizontal size of the target areas (several 
thousands of km) compared with the thickness of the chromosphere under study (a few 
hundred km). The effect of a possible filamentary structure of the atmosphere is however 
much more difficult to assess. 

The only (qualitative) assessment we can make, from a radiative transfer point of view, is 
that photons could perhaps escape more easily from a set of randomly oriented, high-density 
structures embedded in a lower density medium than from a homogeneous atmosphere, 
especially in the case of high density contrast and of a sparse distribution of the filamentary 
structures. Thus, we may expect some modifications in the details of the profiles of the 
more optically thick lines: for instance, a reduction of the depth of self-reversed cores, if 
present. We may also expect that (effectively) optically thin lines should be least affected 
by atmospheric inhomogeneities. 



In the case of more ordered, more closely packed, and thicker structures on the other 
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hand, the outcome could be somewhat different, as shown by iGunar et al.l (120071 ). who 
studied the formation of the hydrogen Lyman lines in prominences, modelled as a set of 
parallel, finite slabs ("threads"), each 1000 km wide: each one, that is, as wide as the entire 
chromosphere we are modelling (Figures 151 and [TTT) . We remark, however, that even in such 
a relatively extreme case, the changes are typically not very dramatic, and observable only 
at high spectral resolution. 

In any case, there is little or no information, to our knowledge, on the effect of a flare 
on the chromospheric fine structure. Thus any more precise assessment on the effect of 
non- uniformity on the emerging profiles is virtually impossible at this stage. 

In Fig. [3] we have shown the sets of optical profiles used in the modeling procedure. 
Since different features of the lines are sensitive to modifications in different parts of the 
chromosphere, the models were constructed following several steps. First, we fixed the 
deepest part of the atmosphere by matching the wings and the Ki minimum of the Ca II K 
line, which are formed in the high photosphere and the temperature minimum region. 

As a second step, we modified the chromosphere and the lowest part of the transition 
region until a satisfactory match was found for the profiles of Ha, Ca II K, and He I 1083 
and 587.6. We recall that the profiles of the first two lines do not depend on the coronal 
back-radiation field. 

The profiles of the two helium lines are formed in two distinct regions: most of the 
radiation is originated in the photosphere, which in the quiet Sun results in a weak absorption 
line at 1083 nm, and no noticeable line at 587 nm. However, already in the active region 
studied in Paper II we found that there is an important chromospheric contribution. This 
contribution becomes, of course, much more important in the flaring atmospheres we are 
studying here, and depends not only on the thermal structure of the high chromosphere and 
the low transition region, between 10 4 and 2.5 x 10 4 K, but also on the coronal EUV incident 
radiation. 

Finally, both ultraviolet lines, the He I line at 58.4 nm and the He II line at 30.4 nm, 
depend on the structure of the low and mid-transition region, from 3 x 10 4 to 5 x 10 4 K for 
the 58.4 line, and up to 1 x 10 5 K for the 30.4 line. As a third step, therefore, we modified 
the structure of this region to obtain a good match for the line fluxes of these UV lines, 
measured with SOHO/CDS. Generally, it was necessary to iterate between steps 2 and 3 
before a final model was obtained. 

We first computed models for regions A and D, marked in Fig. [H using the coronal 
radiation determined in Sec. [3] and the same value of the helium abundance (An e =0.1) 
we used in Paper II. The resulting computed intensity profiles were convolved with the 



-In- 



appropriate instrumental response. 

4.1. Modeling of region A 

For region A, where the observed lines do not change during the considered time interval, 
we obtained the model shown in Fig. [HI which gives a good match between the observed and 
computed lines. In the same figure, we also show, for the line cores, the intensity contribution 
functions, defined as dl/dh = Sexp(— t)x, where / is the intensity, h is the height in the 
atmosphere, S, r and x are > respectively, the source function, the optical depth, and the 
opacity in the line core. 

The observed and computed profiles of the lines formed essentially at chromospheric 
levels are shown in Fig. [9j The bars in the Figure indicate the uncertainty of the observed 
profiles (see also Fig. [3]). It can be seen that the agreement found is very good, well whithin 
the error bars, in all cases. In the case of the optical He I lines, however, we recall that, as 
discussed in Sec. [21 the error bars in the line centers also include the effect of variations of 
the underlying photosphere within the target region, and therefore represent an upper limit 
to the variation induced by chromospheric changes. 

It can also be seen that only the He lines are affected by the EUV radiation, and that 
the computed profiles agree with the observations within the range of variability in the target 
region for both values of the coronal radiation indicated in Table HJ 

We have already mentioned in Sec. [2] that the UV lines (the He I 58.4 and the He II 
30.4) have a large width due essentially to the instrumental profile: therefore the comparison 
with the models can be done only through the radiances, which are shown in Table [2] It 
can be seen that it is possible to obtain line radiances in good agreement with the observed 
ones using the same atmosphere, despite the variation of the coronal radiation of about 50% 
during the considered time interval (see Sec. Sand Tabled]). 

We remark here that the temperature structure in the region of formation of the He II 
30.4 nm line is not constrained by other lines. Thus, for the remainder of the discussion, we 
will not consider a match of the He II Ly-a to be a significant constraint on the atmospheric 
parameters. 
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4.2. Modeling of region D 

For region D, where the chromospheric intensity (Fig. [2]) approximately follows the EUV 
light curves, we were able to find models reasonably matching the observations at t\ and 
t 3 , but it was not possible to obtain a model that would simultaneously reproduce all the 
profiles observed at t 2 . 

In particular, the radiances of the He I 58.4 line always results considerably higher 
than the observed one if we match the profiles of the other lines, and in particular of Ha. 
Conversely, it is possible to find models that reproduce the radiance of the resonance He I 
line, but in this case the other lines, in particular Ho, are poorly fitted. It should be kept in 
mind that it is precisely at this time when the regions of formation of the 58.4 line and of 
Ha overlap most, and therefore the model is best determined (see Sec. I5.ip . 

5. Discussion 

In order to find in region D a model of the flaring atmosphere that matches all the 
observed lines during the considered time interval, and in particular that matches both Ha 
and He I 58.4, we explored the two aspects that most likely affect the helium spectrum: the 
chromospheric helium abundance, An e (Sec. 15.11) . and the spectral distribution and overall 
level of EUV coronal back-illumination (Sec. 15.21) . 

5.1. The effect of helium abundance 

To explore the effect of the helium chromospheric abundance, we computed a series of 
models for kernel D at t 2 , using values of j4H e =0.085, 0.075, 0.065, and 0.055 for T > 6300 K, 
i. e. in the atmospheric layers that strongly contribute to the formation of the considered 
lines. Since a change in the helium abundance affects the hydrostatic equilibrium, models 
with the same T vs. h structure but different abundances give, in principle, different emitted 
profiles not only for the He lines, but for all the other chromospheric lines as well, and for 
Ha and the Ca II lines in particular. Therefore, for each value of Au e we had to compute a 
different model. 

In Fig. [TO] we show the computed radiance of the 58.4 He line and the central intensity of 
the D3, 1083 and Ca II K lines, for the best-matching models for every value of the assumed 
A^ e . The observed values and their uncertainties, as described in Sec. are also shown. 
The variations of Ha around the mean value of its central intensity are of the same order as 
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the analogous variations of the Ca II K line, and thus are not shown in Fig. [TUJ 

The Ca II K and the He I 1083 lines depend very weakly on An e , in the sense that, for 
each value of An e , it was possible to find a model matching the observed profiles within the 
errors. 

On the contrary, the He I D 3 and 58.4 lines strongly depend on A^ e : the 58.4 radiance is 
higher than the observed value (beyond the error bar) for An e = 0.1, and quickly decreases for 
smaller values of A He , crossing the mean observed value at An e ~ 0.07. On the other hand, 
the D3 central intensity increases with decreasing An e , and for A-^ e smaller than ~ 0.075 
significantly exceeds the observed value. In practice, the D3 central intensity approaches the 
photospheric continuum, i. e. the line disappears, for v4 Hc < 0.05. 

In other words, both the He I D 3 and 58.4 nm lines become in fact weaker as An e 
decreases. This result may appear obvious after a superficial analysis neglecting the non- 
linearity of non-LTE problems. In fact, the He I 1083 nm line does not display the same 
behaviour: over the range of abundances explored, it was nearly always possible to find a 
temperature structure fitting both that line and the other chromospheric lines - with the 
exception of the lowest value, 74H e =0.05. 

Considering the error of the measurements, we believe that An c = 0.075 is the best 
compromise that allows reasonably good matches for all the considered lines. We notice that 
a variation of ±14% of A^ e (i- e.: 0.65 < An e < 0.85) still allows a good match between 
observed and computed lines within the error of the observations. 

With this value of ^He, we also determined the models for region D at t\ and £3 UT, which 
are shown in Fig. [TTJ The same figure also shows, for the line cores, the intensity contribution 
functions for the model corresponding to the second maximum of EUV emission, at t 2 . 

A comparison between observations and calculations for the three models is shown in 
Fig. d2l The bars in the Figure indicate the uncertainties of the observed profiles. It can 
be seen that the agreement found is very good, well within the errors, with the same caveat 
regarding the error bars of the optical He I lines discussed in Sec. 14.11 The computed 
radiances for the UV lines are compared with the observation in Table [2j Also in this case, 
the agreement is well within the errors in the observations. 

5.2. The effect of coronal illumination 

In the discussion on the effect of An e on our calculations, we have assumed the values 
of J cor given in Table [TJ Those values, however, have an uncertainty of about 30%, as 
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Region 


Time 


He I 5: 


3.4 nm 


He II 30.4 nm 




(UT) 


Obs. 


Comp. 


Obs. 


Comp. 


A 


16:14 


2.65 


2.77 


16.0 


15.3 




16:33 


2.65 


2.91 


16.0 


14.0 


D 


16:15 


4.34 


4.35 


32.7 


31.5 




16:20 


4.93 


5.29 


24.0 


22.9 




16:35 


3.24 


3.32 


14.6 


15.4 



Table 2: Observed and computed radiances for the He UV lines. Units of 10 4 ergs s 1 cm 2 
sr _1 . 

estimated in Sec. [31 Moreover, the spectral distribution of the photoionizing photons is also 
somewhat uncertain. To assess how these uncertainties affect the profiles of the helium lines, 
we performed a set of calculations for Kernel D at t2, since it corresponds to the second 
emission peak (after the initial impulsive peak), when the differences with the situation 
studied in Paper II are the largest. 

To test the effect of the uncertainty in the overall coronal back-illumination, we com- 
puted the He profiles increasing and decreasing by the ~ 30% estimated uncertainty (Sec. [3]) 
the value listed in Tabled], for the models with An e =0.1, 0.085, 0.075. The results of this 
tests are synthetically shown in Fig. [10] by thin error bars: higher values of / cor produce 
stronger He I lines (deeper optical lines, and higher radiances in the EUV line) . 

For An e =0.1, we found that the computed radiances of the He I 58.4 nm line change 
very little. In fact, a change in J cor much larger than 30% (which, we recall is a "safe" 3-cr 
value) would be required to match this radiance within the error bar of the observed value. 

For AHe=0.075, for which we obtained the best overall fit (Sec. 15.11) we found that 
modifying the value of J cor by 30% does indeed produce significant changes in the calculated 
profiles. It appears, in this case, that the lower value of increases the penetration depth 
of EUV coronal photons, towards regions where the P-R mechanism effectively competes with 
collisional processes in the formation of the helium spectrum, even of the EUV resonance 
He I line. 

To assess to what extent the uncertainty in the spectral distribution of EUV coronal 
photons can affect the results described in Sec. 15.11 we computed the He I profiles using the 
DEM-derived flare spectral distribution described in Sec. 13.11 for the cases with An e =0.1, 
0.085, 0.075. The results are indicated in Fig. [TU] with open dots. 

It can be seen that the changes are small, well within the error of the observations. We 
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can therefore exclude that the discrepancy between calculations and observations in the case 
of A-R e =0.1, noted in Sec. I4.2[ can be due to the uncertainty in the photoionization radiation. 

In order to explore more systematically the effect of the spectral distribution of the 
photoionizing radiation, we also performed tests similar to the ones included in Paper II, 
considering for the spectral distribution different step functions, which are non-zero only in 
a limited spectral range, and zero outside it. The wavelength ranges we considered for this 
tests are: 2.5 - 5.0 nm, 5-10 nm, 10 - 20 nm, 20 - 30 nm, 30 - 40 nm and 40 - 50.4 nm. 

In Fig. [13] we show the intensity in the center of the 1083 line and the radiance of the 
58.4 line, as a function of I cor . We do not show the central intensity of the D3 line, since its 
behavior is very similar to the one of the 1083 nm line. We note that in Paper II we limited 
our tests to EUV radiances 7 cor < 3 x 10 15 photons s -1 cm -2 sr -1 , whereas here we extended 
the tests to 1.2 x 10 16 . 

Several important conclusions can be drawn from this figure. First, the effect is different 
on the two lines: the central intensity of the 1083 line is most sensitive the harder extreme 
of the spectrum (down to the threshold of 2.5 nm), while the reverse is true for the radiance 
of the 58.4 nm line, which is dominated by the longer wavelength bands and it is also nearly 
insensitive to photons below 10 nm. 

For the 1083 and the D 3 lines the relevant parameter, as in the case discussed in Paper II, 
is the ratio between the helium and hydrogen photo-ionization cross sections, which in the 
band 0.5-10 nm is 3 to 5 times larger than at 50 nm. Therefore, coronal photons at larger 
wavelengths are absorbed before they reach the chromospheric heights were these lines are 
formed. Since, as we explained in Paper II, the ratio between the He I and H cross sections 
changes slowly at long wavelengths, the details of the spectral distribution between 20 and 
50.4 nm are of little relevance for the resulting 1083 nm central intensity. 

While these results are qualitatively consistent with those of Paper II, the effect of the 
shortest wavelength photons (in the band 2.5-5.0 nm) is remarkably different. While that 
band had practically no effect on the He I 1083 line computed in the AR model discussed 
in Paper II, we found that, in this flaring atmosphere, these photons may have a role in the 
formation of that line. 

The reason for the difference lies in the higher temperature and densities of the model 
chromosphere of region D: the He I 1083 forms in a region where most neutral or singly- 
ionized metals, whose inner-shell absorption cross sections would otherwise dominate the 
opacity in that band, have already mostly disappeared. 

The He I 58.4 nm line, on the other hand, is formed higher in the transition region, 
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where most hydrogen is ionized, and therefore in this case the relevant parameter is the 
helium photo-ionization cross section alone, which is larger close to the continuum head, at 
50.4 nm. 

Note that, by contrast with the situation studied in Paper II, in this case the much 
larger coronal radiation has some effect on the radiance of the He I UV line, for Au e up to 

0. 1 (Fig. [TUl) . The He II UV line, on the other hand, which forms at higher temperatures, is 
not affected by the coronal radiation, as was the case for the active region prior to the flare 
(see Paper II). 

6. Conclusions 

During a small two-ribbon flare, we obtained simultaneous and cospatial observations, 
including the chromospheric lines Ca II K, Ha and Na I D as well as the He I lines at 587.6 
(D 3 ), 1083 and 58.4 nm and the He II line at 30.4 nm. The EUV irradiance in the ranges 
A < 50 nm and 26 < A < 34 nm has also been measured at the same time. 

We analysed the formation of the He lines after the first impulsive burst of the flare, 

1. e. after 16:13 UT, in two different locations of the flaring area, named A and D. In region 
A the chromospheric radiance is almost constant during the time interval we considered, 
while in region D it approximately follows the EUV light curve with a maximum emission at 
about 16:20 UT (Fig. [2]). A similar behaviour is also detected for the TR and coronal lines 
with temperatures up to 10 6 K. We determined the incident photoionizing EUV radiation 
produced by the corona overlying the target regions, by using information from various 
SOHO instruments (imagers and spectrographs) and we concluded that in region A the 
variability of coronal radiance during the flare is less than a factor 2, while in region D is 
about a factor 5 (Tabled]). 

We determined the spectral distribution of the coronal radiation over the region D at 
16:20 UT, where we observe a peak of total EUV emission, as well as a peak in the hot 
Fe XIX line (Fig. [5]). We found that the largest differences between the pre-flare and the 
flare spectrum occur in the region below 15 nm. 

These differences, however, are not large enough to compensate for the overwhelming 
contribution of longer wavelength photons to the total EUV photoionizing emission in this 
relatively small flare. We therefore simply scaled the spectral distribution of a "typical" 
active region by the appropriate factors listed in the last column of Table [H With this 
coronal radiation, we built semiempirical atmospheric models trying to match our set of 
observables for regions A and D, using the same value of the helium abundance (An e =0.1) 
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we used in Paper II. The He lines and continua are self-consistently computed in the radiative 
transfer calculations. 

For region A we were successful in obtaining line spectral intensities and radiances in 
good agreement, within the errors, with the observed ones at the considered times, using a 
constant atmosphere and the coronal radiations determined in Sec. [3j Taking into account 
that the formation of the helium lines in this case still falls in a regime where the P-R 
mechanism is important (Fig. [H]), this result makes us confident that our assumptions for the 
coronal radiation and its spectral distribution are acceptable. 

For region D, we were able to find models producing profiles reasonably matching the 
observations at 16:15 and 16:35, but it was not possible to obtain a model capable of simul- 
taneously reproduce all the profiles observed at 16:20, when the flare emission is maximum. 
In particular, it was not possible to simultaneously match Ha and the He I 58.4 line. We 
therefore explored the two aspects that most likely can have an effect on our calculations: the 
chromospheric helium abundance, An e (Sec. 15.11) . and the spectral distribution and overall 
level of EUV coronal back-illumination (Sec. 15.21) . 

Since a change in the helium abundance affects the hydrostatic equilibrium, for each 
considered value of A^e we had to compute a different model. We used values of AH e =0.085, 
0.075, 0.065 and 0.055 for T > 6300 K and we showed that the He I D 3 and 58.4 lines 
strongly depend on An e , both becoming weaker as A^ e decreases. In particular the D 3 line 
disappears for An e < 0.05, setting a lower limit for the A He value. On the contrary, over 
the range of abundances explored, it was always possible to find an atmospheric structure 
matching He I 1083 and the other chromospheric lines. Taking into account the error of the 
measured lines, we regard A^ = 0.075 as the best compromise that allows reasonably good 
matches for all the considered lines and we found that a variation of ±0.010 (or: ±14%) of 

is still acceptable within the error of the observations, with the average value of coronal 
illumination determined in Sec. 15.^1 

On the other hand, an analysis of the effect of the coronal illumination at a given value of 
Aft e (in Sec. I5.2l we considered A He =0.1, 0.085 and 0.075) led us to exclude that uncertainties 
in the determination of this parameter could be responsible for the discrepancies at v4H e =0.1. 

We also determined new models for region D at 16:15 and 16:35 UT using An e =0.075. 
Thus, for these two times we have models matching all the observations for both A^ e = 0.075 
and A^ e = 0.1. In these cases, as in the case of Region A (Sec. 14. lj) the incomplete overlap 
between the formation regions of Ha and He I 58.4 (see for instance Fig. [HD does not allow 
a determination of An e with an uncertainty better than a factor of two, as in Paper II. 

In summary, our analysis of the observed chromospheric spectra, is consistent with 
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AHe=0.075; furthermore, we were able to estimate the uncertainty range as 0.065-0.085 in 
at least one favor able case. The value A^ P =0.085 - i. e. the "canonical" photospheric helium 



abundance (e. g.: lAsplund et al.ll2005l )) - could therefore be marginally consistent with our 
data, since it results in He features which agree with the observations at the limit of our 
estimated errors. 

A final note on measuring the helium abundance in the solar chromosphere: The set 
of lines we have used, together with an estimate of the EUV coronal back-radiation, is the 
minimum necessary to estimate An e in chromosphere. No meaningful estimates, for instance, 
can be obtained using only the optical He I lines or the EUV resonance lines alone. But (as 
in Paper II, or in all cases considered here, except region D at 16:20 UT), even this set of 
lines may not be sufficient for an accurate determination of An e . In particular, while the 
optical He I lines are usually well constrained by typical chromospheric lines such as Ca II K 
or Ha, the resonance He I 58.4 line must also be constrained by lines forming in the lower 
transition region. In the high pressure environment of region D at the secondary peak of the 
EUV flare emission, this role is played by Ha. In other circumstances, additional constraints 
provided by other lines would be needed: intensities from other upper-chromospheric or 
lower-TR lines (e. g.: HI Ly-a); intensities from other He I and He II lines (e. g.: He II 
25.6 nm); spectrally resolved profiles of the He I 58.4 and He II 30.4 lines. With the advent 
of more recent space missions, such as Hinode and its EUV spectrometer (EIS), and with the 
continued availability of the SOHO spectrometers (CDS and SUMER), there is now in fact a 
window of opportunity for obtaining stronger constraints for the kind of modeling described 
in this work. 
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Fig. 2. — Light curves in the EUV and soft-X bands (1-minute resolution), before and during 
the C-class flare studied in this paper, compared with optical measurements. Top panel: 
EUV Irradiance at 1 AU, from SOHO/SEM (black: A < 50 nm; grey: 26 < A < 34 nm); 
Middle panel: Soft-X irradiance from GOES-10 (0.05-0.4 nm band); Bottom panel: Ha 
intensity observed at the Dunn Solar Telescope (DST) in the two representative regions 
analyzed in this paper (see Fig. [1]). The vertical arrows mark the times of the full-disk 
images in the 19.5 nm band taken by SOHO/EIT during the CME watch observing program, 
while the vertical dashed and dotted lines mark the times of the spectra of regions A and D, 
respectively. 
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Fig. 3. — Profiles for the lines Ha, Ca II K, He I 1083.0 and 587.6. Top panels: Observations 
in a quiescent, reference region (dashed lines) and in region A at 16:14 UT (solid lines). 
Bottom panels: Observations in region D at 16:15 (dotted lines), 16:20 (solid lines), and 
16:35 UT (dashed lines). Error bands for the profiles in region A at 16:14 UT and region D 
at 16:20 UT are also shown. Note: telluric lines have been removed from the He I 587.6 
profiles, except in the case of the reference region. 
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Fig. 4. — Sequences of CDS/NIS residual spectroheliograms showing, for about 40 minutes, 
the evolution of the event from its onset (at about 16:01 UT). The images obtained from the 
profiles of the He I 58.4 nm are very similar to those of the He II 30.4 nm line, and thus are 
not shown here. Circles mark the two regions of interest (A and D: left-most and right-most 
circles, respectively). 
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Fig. 5. — Top panel: Light curves of residual irradiances in CDS lines; the EIT enhancement 
over the flaring region is also shown. Bottom panel: Light curve of excess SEM irradiance 
during the flare (black, thick line, and empty circles), compared with the fit (grey line, and 
filled circles) obtained by combining the CDS and EIT light curves. The vertical dashed and 
dotted lines mark the times of the spectra taken in regions A and D, as in Fig. [2l 
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Fig. 6. — Sequences of images showing the variation of the enhancement of EUV illumination 
in the A < 50 nm band, in the flaring area. Radiance maps shown with the same format as 
in Fig. HI 




log T [K] Wavelength [nm] 

Fig. 7. — Left: DEM distributions for the "reference" AR atmosphere (black line) and 
the flaring region D at 16:20 UT (grey line, yellow in the on-line version). The observed 
intensities (multiplied by the factor DEM(T cff )// ca i c ) are also shown for the various CDS 
lines: He I 58.4 nm (triangles), He II 30.4 nm (x symbols), O V 62.9 nm (diamonds), 
Fe XVI 36.1 nm (squares), and the blend at 59.2 nm (filled circles). Right: Spectral distri- 
butions from the above DEMs, excluding all He I and He II lines and continua, compared 
with the "standard" AR emission times the multiplicative factor of ~ 7.7 given in Table [1] 
(dashed histogram). 
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Fig. 8. — Region A: temperature versus height distribution (solid line) of the atmospheric 
model obtained for 74H e =0.1. The model for the active region studied in Paper II (dotted 
line) is displayed for comparison. Also shown are the intensity contribution functions, in 
arbitrary units (solid, thick lines), of the different spectral features we used to build the 
model. The He I 587.6 contribution function is proportional to the one of the infrared He I 
line. 
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Fig. 9. — Region A: comparison between observed (thick grey line) and computed profiles 
(black lines). The solid and dashed lines indicate the profiles computed for the radiation 
determined at 16:14 and 16:33 UT, respectively. 
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Fig. 10. — Computed radiance (10 4 ergs s _1 cm -2 sr _1 ) and central line intensity (10 7 ergs 
s _1 cm -2 sr _1 nm" 1 ) versus ^He- Radiance of the 58.4 line and the line center intensity of the 
D 3 , 1083 and Ca II K lines are computed for the best-matching models obtained for region 
D at t 2 , with I cor given in Tabled] and the scaled spectral distribution of the active region 
(crosses). Open dots indicate values obtained with the flare spectral distribution described 
Sec. 13.11 The thick line indicates the observed value and the dotted line its estimated 
uncertainty. The thin line indicates the continuum level for D3. 
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Fig. 11. — Region D: temperature versus height distribution of the atmospheric models 
obtained for v4 Hc =0.075 at the considered times (dotted line: 16:15 UT, solid line: 16:20 UT, 
dashed line: 16:35 UT). Central intensity contribution functions of various lines are also 
shown, as in Fig. El 
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Fig. 12. — Region D: comparison between observed (thick grey line) and computed profiles 
(black line). The profiles at 16:35 and 16:20 UT are offset by an arbitrary, constant amount 
along the ordinate. 
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Fig. 13. — Radiance for the 58.4 nm line (top), and central intensity of the 1083 line (bottom), 
as a function of J cor , for different step- like distributions of the field. Same units as in Fig. [TD] 
for the ordinates; I cor given in 10 14 ph. s _1 cm -2 sr -1 . The vertical line marks the value 
of I cor at £2, while the horizontal lines indicate the mean value of the measured 58.4 nm 
radiance (top panel) and 1083 nm central intensity (bottom panel). 



